Exact RG computation of the optical conductivity of graphene 
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The optical conductivity of a system of electrons on the honeycomb lattice interacting through 
an electromagnetic field is computed by truncated exact Renormalization Group (RG) methods. 
We find that the conductivity has the universal value tt/2 times the conductivity quantum up to 
negligible corrections vanishing as a power law in the limit of low frequencies. 

PACS numbers: 71.10.Fd, 72.80.Vp, 05.10.Cc, 05.30.Fk 



Among the remarkable properties of graphene 1 , the 
optical conductivity is of special interest. Recent 
experiments 2 found that the conductivity in monolayer 
graphene is essentially constant in a wide range of fre- 
quencies between the temperature and the band-width. 
The observed value of the conductivity is equal, within 
experimental errors, to oo = 7re 2 /(2/i), a universal value 
that only depends on fundamental constants and not on 
the material parameters like the Fermi velocity. This 
fact can be nicely explained in terms of the standard 
graphene's model of massless Dirac particles in 2+1 di- 
mensions: in this case, by neglecting interactions and 
disorder, Ref. 3 predicted a value of the conductivity 
073(0;) at inverse temperature ft and frequency uj satisfy- 
ing lim^-^o lini/3-i.oo o~p{uj) — <7q. The inclusion of lattice 
effects does not change the value of this limit 4 . 

This remarkable agreement with a theoretical value 
computed by neglecting many body interaction is, how- 
ever, surprising and needs an explanation 5 . Indeed, the 
strength of the interactions in graphene is measured by 
the ratio a = ^ ~ 2.2 (e is the electric charge and 
vo the Fermi velocity), which is 300 times larger than 
the usual fine structure constant. The effects of the in- 
teractions are clearly seen in experiments on the Fermi 
velocity 6 . Therefore, why is not there an essential many 
body renormalization of the optical conductivity, too? 

On the theoretical side, a theorem proved in Ref. 7 es- 
tablishes that the conductivity of electrons hopping on 
the honeycomb lattice and interacting via a weak Hub- 
bard interaction is equal to o~o in the limit uj — > 0. Note 
that, even if dimensionally irrelevant, the interaction can 
produce finite many body renormalizations: for instance, 
the Fermi velocity is renormalized by the interaction. 
Therefore, the universality of the conductivity is a non 
trivial statement, following from an exact cancellation of 
all the many body corrections. 

It is, however, believed that the interaction in clean 
suspended graphene is not at all short-ranged as the Hub- 
bard interaction (no screening), so that a more realistic 
description of the clean system requires the inclusion of 
the long-ranged electromagnetic (e.m.) interactions. In 
the case of static Coulomb interactions, Ref. 8 predicted a 
logarithmic renormalization of the Fermi velocity, namely 
v(q) = uo(l+ § log |) where q is the momentum mea- 
sured from the Fermi points and e is the bandwidth. 



First attempts to include the effects of a Coulomb po- 
tential on the conductivity 9 led to the conclusion that 
the interaction radically changes its behavior, that is 
lim w _>.o o~{uS) — 0, where o-(lu) is the conductivity in the 
limit of zero temperature. Later, Ref. 10 ' 11 obtained the 
qualitatively different result lim w _i.o <T ( W ) — Co, based on 
scaling arguments. In particular, Ref. 10 found the for- 
mula 



tr(u)=<r [l + O( 



1 



log(e/w) ' 



(1) 



Note that the inverse logarithmic correction in Eq.(l) is 
a consequence of the logarithmic divergence of the Fermi 
velocity, and should be read as 0{q.vq/v{uj)). As pointed 
out in Ref. 5 , this correction is in general larger than the 
experimental error 2 . Ref. 12,13 proposed that the way 
out from this apparent contradiction should be found 
in the constant in front of the inverse log corrections, 
whose correct value should be much smaller than the one 
computed in Ref. 10 . However, Ref. 14 raised objections 
against the new value proposed in Ref. 12,13 , because the 
regularizations used in these works can produce unphys- 
ical results. The disagreement between the big (inverse 
log) corrections to the conductivity and the experimental 
data suggested 5 to phenomenologically postulate a Fermi 
liquid description of the interacting system: this assump- 
tion implies that the universal conductivity is reached at 
low frequencies polynomially fast (i.e., as ~ uj 2 ) but is in 
contrast with the experiments in Ref. 6 . 

Eq.(l) was derived by assuming that the electrons in- 
teract via a static Coulomb interaction: however, the log- 
arithmic increase of the Fermi velocity shows that the as- 
sumption of instantaneous interactions becomes unphys- 
ical at very low energy scales 15 . Therefore, the use of 
Eq.(l) and of the divergence of the Fermi velocity to pre- 
dict the universality of the conductivity as uj — > is ques- 
tionable. The unbounded increase of the Fermi velocity 
is absent in the case that the interaction with the e.m. 
field is introduced via the Peierls substitution in order 
to preserve gauge invariance. It is well known 16 ' 17 that 
in this case the Fermi velocity stops flowing at the speed 
of light c and Lorentz symmetry spontaneously emerges 
in the infrared. We compute the optical conductivity at 
imaginary frequency uj in a lattice gauge invariant model 
for graphene using truncated exact RG methods. We find 
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an expression that is qualitatively different from Eq.(l), 
namely 18 



a( W )=<7 [l + O(-log-)] 



(2) 



which is very close to the universal one at low frequencies, 
up to a really negligible power law correction, compati- 
ble with the experimental results in Ref. 2 . The wlogw 
dependence of the correction is not necessarily optimal, 
it may just be a byproduct of our estimates. 

We derived Eq.(2) under the assumption that the val- 
ues of the bare parameters are sufficiently close to the 
infrared fixed point (i.e., the bare Fermi velocity vq is 
sufficiently large). The extension of its validity to real 
frequencies and to a larger range of parameters, includ- 
ing those measured in actual graphene's samples, requires 
a microscopic justification that is quite difficult in view 
of the strength of interactions in graphene 19 ; of course, 
this is a caveat that applies to all the approaches based 
on expansions, resummations and truncations. In any 
case, it is reassuring to see that encoding a fundamental 
physical principle like gauge invariance into the model 
is sufficient to obtain results that are in good qualitative 
agreement with the experimental data, in particular with 
the observed dramatic increase of the Fermi velocity and 
with the universality of the conductivity up to negligible 
power law corrections at low frequencies. 

The model we consider was defined in detail in Ref. 17 . 
Let us just remind here the main definitions. The grand- 
canonical Hamiltonian at half-filling is H = Hq + He + 
Ha, with 

H = -t T Y a+ b- . e ^So 1 ^As+s5 j fi) ds + cx _ 

t-~< Z-^i x,<r x+Sj,a 

xeA A o-=T4- 
.7=1,2,3 

the gauge invariant nearest neighbor hopping term (here 
t is the hopping strength, 5j the nearest neighbor vectors 
and a ± ,b ± the creation/annihilation operators of elec- 
trons sitting at the sites of the A or B sublattice of the 
honeycomb lattice), 



S,yeA A UA B 



where e is the electric charge, ipp is a regularized version 
of the static Coulomb potential and ng the electron num- 
ber at site x. Finally, Ha is the energy (in the presence 
of an ultraviolet cutoff) of the three-dimensional photon 
field A = (A, A 3 ) in the Coulomb gauge. Units are fixed 
in such a way that the speed of light c = 1. Note that 
the interaction with the quantum e.m. field is introduced 
via the Peierls substitution in order to preserve Gauge in- 
variance. 

Proceeding as in Ref. 7 , where we computed the conduc- 
tivity in the case of short range interactions, we define a 
"space-time" three-components vector Jp l([i , n — 0,1,2, 



with 



SeA A 



SeA B 

<r=U 



the density operator and Jp,i,Jp t 2 the two components 
of the paramagnetic current 



L=iety e-^lrfJat bZ t 



xGA 



x+8j : a 



x+dj 



where r] J ~ = (1 — e vp j ) / (ipb~ j) . Let also p = (u!,p), with 

uj the Matsubara frequency, and K^fjp) is the current- 
current response function, i.e., the Fourier transform of 
lini/3-s-oo (J Xj([t ; J y ,v)o- 

We are interested in the conductivity, defined via Kubo 
formula as 4,7 (here l,m — 1,2): 



o-lm{u) = - 



2 1 



K lm (uj,d) + A lm (6) 



where 3^3/2 is the area of the hexagonal cell of the hon- 
eycomb lattice and 



A im (p) 



lim —r y 

x£A 
.3=1,2,3 



with a sj = -#tY,*(4j s+ i j , a + K+s il „ a 2J thc dm ~ 

magnetic tensor. 

The current-current response function can be com- 
puted via the generating functional that, in the Feynman 
gauge, reads 



W h .(J,A) 



= J P(#) J 



P h *{dA)e 



V(A+J,i/,) + M>,A) 



(4) 



In 



which has been studied in great detail in Ref. 
Eq.(4): (i) tp^ are Grassman spinors (of thc form 
tp = (a,b), with a and b the electron fields associated 
to the two sublattices of the honeycomb net) and P(dip) 
is the fermionic gaussian integration with propagator 



ff(k) = - 



j_ / ik v n*(k) 

Zq \ VoQ(k) ik a 



(5) 



where Zo = 1 is the free wave function renormaliza- 
tion, v = |i is the free Fermi velocity and fl(k) = 

| 2 3 e lk ( Sj ~ Sl ^ is the complex dispersion relation. 

Note that g(k) is singular only at the Fermi points p r F = 
(0, 2 f", r ^|), where r = ± is the valley index. More- 
over, A fl (p), n = 0, 1, 2, are the Fourier transform of real 
gaussian variables and Pf {dA) is the gaussian integra- 
tion with propagator «v(p) = 5 Ail/ (2|p|)^ 1 X[ft.«,o] 
where X[h*,o] is a smooth compact support function that 
acts both as an ultraviolet cutoff on scale |p| ~ 1 and as 
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an infrared cutoff on scale |p| ~ 2 h (to be eventually re- 
moved) . Finally V is the interaction whose explicit form 
can be easily inherited from H 17 . The current-current 
response function can be obtained by taking the limit 
h* — > — oo and by deriving twice with respect to the ex- 
ternal field J and then setting J = A = 0. Field- field 
correlations or field-current correlations can be obtained 
similarly, by suitably deriving with respect to the exter- 
nal fields A and/or J. Note that in writing the generat- 
ing functional as in Eq. (4) we exploited gauge invariance 
and, more precisely, the equivalence between the Feyn- 
man and the Coulomb gauges. Another crucial conse- 
quence of gauge invariance is the following equation 



= 



d 



J + da,\e iea ) 



a=0 



(6) 



By performing derivatives with respect to the external 
fields, this equation also implies a sequence of exact lat- 
tice Ward Identities, valid for each finite choice of the 
cutoff scale h*. In particular, proceeding as in Ref. 7 and 
defining p° = —ioj, the current-current response function 
satisfies the Ward Identities Eu=o?'''^o(p) = and, 
for j = 1,2, 



^K^ j (p) = -J2p i A ij (p). 



(7) 



An immediate consequence of Eq.(7) and of the continu- 
ity of K^, u (p) in p = (proved at all orders of renormal- 
ized perturbation theory 17 ) is that, for i,j e {1,2}, 



(Jij (lj) 



2 1 



[Xy(w,0)-.Ky(0,i 



(8) 



see 7 for the simple argument leading to Eq.(8). 

The generating function (4) can be computed by ex- 
act RG methods 17 , which allowed us to prove that the 
response functions can be written in terms of a renor- 
malized perturbation theory that is finite at all orders 
in the effective coupling constants, with explicit bounds 
on the n-th order contributions. In particular, after the 
integration of the degrees of freedom corresponding to 
momenta larger than 2 h , h < 0, we rewrite (setting for 
simplicity A = 0): e w "*W°) = 

= j J] P(d^^)P(dA^' l ))e v(h, (^ < - h, ^ < - h, + J ) , 

' r=± 

where P(dipi- h ^) and P(dA(- h ^) have propagators 



# ft) (k') = - 



Xh(k') 



ik 



Z h \v h n(k' + F F ) 



v h n*{k' + p F ) 

ik 



(9) 

and wlf l) (p) = S^(2\p\)- 1 X[h-,h](p), where: (i) Xh(k') 
is a smooth cutoff function vanishing for momenta larger 

than \k'\ ~ 2 h ; (ii) X[h*,h](p) = Xh(p)-Xh'(p)\ Z h ,v h 
are the effective wave function renormalization and Fermi 



velocity at scale h. Moreover is the effective poten- 
tial, expressed by a sum of monomials in ip(- h \ A^- h ^ of 
any degree: 



V {h \^Z~ h ^ h \A^) 



dp 



[7W p5 (</>) A(<h) 



(2tt) 3 1 j p>p m.p 

2V,,,i^ p i^ ) ] +m^>(v/z^^U^)) 



(10) 



where KV^ (y r Z h ~ip ( -- h) , A^- h) ) is the irrelevant part of 
the effective potential (sum of all the terms with more 
than three fields) and 



*(<h) - 



(3L 2 



r,cr,k' 



with r°(fc) = 1 and 



— i e ik ( S 3~ S l) 




We can summarize the previous discussion by saying that 
after the integration of the degrees of freedom corre- 
sponding to momenta > 2 h , we get an effective theory 
that is qualitatively very similar to the original one, mod- 
ulo the renormalization of a finite number of effective pa- 
rameters, namely the Fermi velocity Vh, the wave func- 
tion renormalization Zh, the vertex function Z,^ and the 
photon mass v^^. The discrete rotational symmetries of 
the model imply that Z^ = Z^ and = i>2.h- 

These parameters verify suitable flow equations well 
defined at all orders in the renormalized expansion: this 
is an instance of the renormalizability at all orders of 
the model. Note that the effective charges at scale h 
are: e {) . h = eZ^/Z h and e^ h = eZ ( h l) / (Z h v h ), for i e 
{1,2}. Thanks to the Wis induced by Eq.(6), we proved 
in Ref. 17 that v^h = 0(e 2 2 h ) (i.e., the effective photon 
mass vanishes in the infrared) and that the beta function 
for the effective charge is asymptotically vanishing, i.e., 



eo,h 



7 (0) 

J z h " 

; 

' Z h v h 



e{l + 0(el h )) , 
= e(l + 0(e^)) 



(12) 



and eo.-co = ei,-oo = 62.-00- Moreover, see 17 , the wave 
function renormalization diverges in the infrared, while 
the effective Fermi velocity increases up to the speed of 
light, both approaching their limits with an anomalous 
power law: 



-■qh 



1 - V h 



2 fjh 



2ei 



5tt 2 



(13) 

+ Oiei^) the 



with n = + Oietvo) and r\ 
two critical exponents. 

The above integration procedure leads to an expan- 
sion of the conductivity in terms of powers of e^.h', such 
renormalized expansion is a resummation of the naive 
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perturbative expansion in e. It must be stressed that 
there is a big difference between these two expansions: 
while the one in e^.h is order by order finite (with explicit 
bounds on the growth of the n-th order contributions 17 ), 
the naive one in e is plagued by O (\og n uj) divergences 
at order n. Therefore, the truncation of the renormal- 
ized expansion is expected to give much more accurate 
predictions than the naive one. 

By truncating the exact RG expression for the conduc- 
tivity at one loop, we get contributions from the bubble 
diagrams in Fig.l, 



h 




h' 

FIG. 1. The one-loop bubble diagram contributing to the 
longitudinal conductivity au(uj). The labels h, h indicate 
the scale of the two loop propagators, which depend on the 
dressed Fermi velocities v h , v h i and on the effective wave func- 
tion renormalizations Z h ,Z h i. A summation over h, b! is un- 
derstood. The big dots correspond to dressed vertex functions 
with h — max{/i, h'}. The external momentum flowing 
in the wavy lines is p = (w, 0). 

which give (defining cr = § x) : 

16 1 ^ r dk f dk> {zff 
a ° " 3V3^^ i 2tt J B \B\ Z h Z h , ■ 1 J 

■^{ri(k')gi h \k')ri(k^ r h \^ + (w,6)) - . 9 W(k')] } • 

where: (i) B is the first Brillouin zone and \B\ = 
8tt 2 /(3%/3) its area; (ii) h = max{/i, ft'}; (iii) T l r (k') = 
T % (k' + pp.); (iv) g^ is the effective propagator on 
scale h, given by the same expression as Eq.(9) with 
Xh(k') replaced by the smooth compact support function 
fh(k') ■= Xh(k') - Xh-i{k'), which is non vanishing only 
if k 1 is on scale h, i. e., 2 h - x < \k'\ < 2 h+1 . Note that 
the effective parameters Z^\z h , vh entering Eq.(14) are 
all functions of e: if we expanded them in e we would 
recover infinitely many graphs of the naive perturbation 
theory, all plagued by logarithmic divergences. Note also 
that Eq.(14) is not simply the "bubble graph" with the 
dressed propagator and vertices: e.g., if one thinks of the 
dressed propagator with momentum k as being obtained 
by resummations of the chain of self-energies, one has to 
take into account that the scales of the momenta flow- 
ing inside such self-energy sub-diagrams are higher than 
the scale of k, according to the rules of exact RG (which 



avoid the problem of overlapping divergences and, corre- 
spondingly, the emergence of nl factors at higher orders). 

The computation of Eq.(14) can be explicitly per- 
formed, by making use of Eq.(12) and by carefully ex- 
ploiting symmetry cancellations that make the apparent 
logarithmic divergence of Eq.(14) finite, see Appendix A 
for details. The result is Eq.(2). 

Our analysis is based on a truncation of the exact RG 
equations, and the question of how to generalize it to the 
full RG expansion is a very interesting and important 
theoretical problem; so far, we succeeded in performing 
the full RG computation only in the case of short range 
interactions 7 . Another important open problem is to un- 
derstand the analytic extension of the conductivity to 
real frequencies. 

In conclusion, we considered a model of electrons 
on the honeycomb lattice interacting via a quantized 
photon field, previously investigated in Ref. 17 . The 
coupling with the e.m. field is introduced via the Peierls 
substitution in order to preserve gauge invariance. We 
showed that at low frequencies the conductivity is equal 
to the universal value <r up to corrections O(wlogw), 
which are much smaller than the 1 / log u> corrections 
found for static Coulomb interactions. Our results 
are a priori valid close to the infrared fixed point and 
the extension of their validity to a larger range of 
bare parameters (including those measured in actual 
graphene's samples) is based on a phenomenological 
assumption. Still, it is reassuring to see that it is enough 
to encode gauge invariance in a microscopic model for 
clean graphene to recover good qualitative agreement of 
the predictions with the experimental data. 

Acknowledgements We acknowledge financial sup- 
port from the ERC Starting Grant CoMBoS-239694. 

Appendix A: Computation of the conductivity 

Using the definition of gi h ^ into Eq.(14), we can rewrite 
(shifting the momenta by p r F and neglecting for consis- 
tency terms of order 0(e 4 ), which should be combined 
with the two-loops contributions) 

a u (uj) _ 4 1^ rdkp f d fj ( Z h) 2 fh(k')f h '(k') 

*o "~ tt 2 w h j^J 2^ J B z h z h , e o + v i\n r (kT ' 

'-k (k + i^)K r (fc')| 2 + v h v h ,Re(n 2 r (k')al r (k')) 

(k + w) 2 + vl,\n r {k')\ 2 

| fc 2 |a t , r (fc')| 2 -^^Re(^(fcOa 2 r (fcQ) l 
k 2 +vl\Q r (k')\ 2 \ ' 

where Q r (k') = Cl(k' + Pp) and 



■5 



Note that fi r (fc') = ik[ + rk' 2 + 0(|fc'| 2 ), a hr {k') = 
i + 0(\k'\) and a 2 , r $) = ~ r + °(W\)- The "relativis- 
tic approximation" consists in replacing O r (fc'), a\ tr (k') 
and a2, r (k') in Eq.(Al) by ik[ + rk' 2 , by i and by — r, re- 
spectively. By performing this replacement, it becomes 
apparent that the r.h.s. of Eq.(Al) behaves dimension- 
ally as i/ rf3fc/ [ (Fqb? - (Fp"]' wnicn is logarithmically 
divergent as oj — > 0. In order to prove the finiteness of 
cr a (oj) in the low frequency limit, it is necessary to ex- 
ploit cancellations, which follow from the use of a Ward 
Identity combined with an essentially explicit computa- 
tion (using the residues' theorem to integrate ko out) of 
the r.h.s. of Eq.(Al). 

Let e be a small but finite fraction of the bandwidth 
(say, e = t/10) and let us distinguish the contributions to 
the integral coming from the region t>o|f2 r (fe')| > e from 
those Wo|O r (fc')| < e. The former correspond to non- 
singular contributions, which can be estimated as fol- 
lows: we expand in Taylor series the expression in square 
brackets up to O(oj 2 ); we note that the term linear in to 
is vanishing by parity in ko; we bound dimensionally the 
term quadratic in w as: 

1 f v 2 oj 2 OJ 

(const.)- J dk jjk'k' (k2+ V ^ {k/)2)2 < (const.)^ . 

(A2) 

Let us now look at the the terms coming from the re- 
gion wo|O r (fc')| < e. Note that the contributions from 
r = + or r = — are equal among each other, thanks 
to the symmetry under valley exchange. We intro- 
duce some shorthands: we define Ah = u/j|f2 + (fc')| and 
Wh,h' = VhVw'Re{p? + {k')a1_ + {k')) , so that, by performing 
the integral over ko using the residues' theorem, we can 
write the contribution to the conductivity coming from 
the region vo\Q r (k')\ < e as 



dk 



«>>2 



^ W J\k>\<e/v7" h y <0 Z hZh> 



■h(k')f h ,(k>) 



+ 



+ 



A h {A h - ioj)\ ai . + (k')\ 2 + W h , h > 
2A h [oJ 2 + 2iojA h -(A 2 -A 2 ,)] 

A h ,(A h , + iw)\ai,+ (k')\ 2 + W h ,h> 
2A, l ,p-^A,, + (A 2 -A 2 ,)] 
A 2 K + (fc')| 2 + W hM A 2 ,|« t . + (fc')| 2 + \Y h ,,- 
2A,(A 2 -A 2 ,) 



(A3) 



2A^(A 2 



+- w hih > \ 
A 2 ) J 



Let us first consider in Eq.(A3) the terms with the inte- 
grand proportional to Wh.h', which can be rewritten as 
(defining V h . h < = A h + A h >) 



I 



|fc'|<e/u 



dk^_^ 

h,h'<0 



(zf) 2 W hJl , -4ojf h (k')fh>(k') 
Z h Z h , A h A h , V h , h >(uj 2 + V 2 h ,) ' 



(A4) 



Now note that W h ,h> = v h v h ,[{k' 2 ) 2 - (k[) 2 + 0(|A?'| 3 )]: 
therefore, the term proportional to (k' 2 ) 2 — {k[) 2 is zero 
by symmetry, and we are left with a contribution dimen- 
sionally bounded as (using that A/,, Ah> and Vh,h' behave 
dimensionally as ~ k' close to the singularity) 



(const.) oj 



rE/V ° V (k') 3 . , OJ , £ 



L " k ' k 'w 



- , 2nlv2\ - (const.) — log ■ 
(A5) 



for oj -C £. We are left with the terms obtained by re- 
placing WhM with in Eq.(A3), which are given by 

4 r (Z^) 2 
- w /- & E -^+Mk')h>{k')- 

77 JW\<e/v h j^ Q Z h Z h , 
K+(k')\ 2 

\n+(k')\(v h + v h ,)[uj 2 + ( Vh + v h ,) 2 \n + (k')\ 2 ] ' 

Using the Ward Identity Eq.(12) and the rewritings 
fi+(#) = ik[ +k' 2 + 0(|fc'| 2 ), a h+ (k') = i + 0{\k'\) 
and d2,+ (k') = —l + 0(\k'\) (valid close to the singu- 
larity k! = 0), this is equal (up to terms coming from the 
"non-relativistic" parts of Cl + (k') and ai i+ (k'), which are 
bounded as in Eq.(A5)) to 



-E / 



i/«o Z- 

dk' h 



fh(k')f h '(k') 



z h z h , (^ ) [ w2 + (^ )24 (fc0 2 ] 

(A6) 

By the compact support properties of ft, the integrand 
is non vanishing only of \h — h'\ < 1, in which case 

^^ = l + 0(e 2 ), -^^l + (e 2 ), 
uniformly in h. Therefore, Eq.(A6) is equal to 



f — 

Jo * 



Aoj 2 2e. , oj 2 

P + 4(fc0 2 ] = ^ arCtan( - } = + } ' 

(A7) 

up to terms bounded by (1 —v )oj/e and terms bounded 
uniformly in oj by 



(const. )e 2 / dk'—— 
Jo ( fe 



) 2 +^ 



< (const.)e 2 . (A8) 



Of course, these correction terms should be neglected, for 
consistency, because Eq.(14) is obtained by truncating 
the exact RG expansion at one loop. Putting all together 
we get Eq.(2). 
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